function x=TikhonovHilbert(alpha, n)
H = hilb(n+1);
s = zeros(n+1, 1);
for j = 0 : n
    s(j+1) = (-1)^j;
end
t = zeros(n+1, 1);
for j = 0 : n
    for i = 1 : j
        t(j+1) = t(j+1) + (-1)^(i+j) / i;
    end
end
r = s * 0.69315 + t;
B = H' * H + alpha .* eye(n+1);
y = H' * r;
x = gausslineqn(B, y, n+1);

end